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Abstract 

Diffusion of Cu adatoms and dimers on Cu(lll) and Ag(lll) surfaces is analyzed 
basing on ab-initio surface potentials. Single adatom diffusion is compared with 
dimer diffusion on both surfaces. Surface geometry makes the adatoms jump alter¬ 
nately between two states in the same way in both systems, whereas dimers undergo 
more complex diffusion process that combines translational and rotational motion. 
Small difference in the surface lattice constant between Cu and Ag crystals results in 
a completely different energy landscape for dimer jumps. As an effect the character 
of diffusion process changes. Homogeneous Cu dimer diffusion is more difficult and 
dimers rather rotate within single surface cell, whereas diffusion over Ag surface is 
faster and happens more smoothly. Temperature dependence of diffusion coefficient 
and its parameters: energy barrier and prefactor is calculated and compared for 
both surfaces. 


1 Introduction 


Surface diffusion coefficient is one of the factors that decide about the crys¬ 
tal growth process. In order to understand formation of nanostructures it is 
necessary to understand the underlying microscopic processes, such as adsorp¬ 
tion, desorption and diffusion. A lot of studies focus on dynamics of a single 
adatom mm- However, diffusion of small clusters, e.g., dimers and trimers, 
which occurs with lower barriers gElEl , can be also signihcant and can af¬ 
fect island formation.In particular on metallic surfaces even diffusion of larger 
clusters can be very effective P7] . 
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Diffusive motion of atoms and clusters at the surface is studied experimentally 
by scanning tunneling microscopy (STM) and field ion microscopy (FIM), 
both used for direct imaging of individual adatoms piHlllUllllj . With the help 
of other techniques, such as held emission microscopy (FEM), photoemission 
electron microscopy diffusion coefficient can be measured. 

As for theoretical approaches, the hrst step is always calculation of the adia¬ 
batic potential, that means all the adsorption sites and energy barriers between 
them [T^13fl4lll5fl6] . Then these data are used to analyze the diffusive mo¬ 
tion. It can be done by means of simulation methods, such as kinetic Monte 
Carlo (kMC) or molecular dynamics |15lll6lll7lll8lll9ll2nj . The latter solves the 
equations of motion for adatoms in a given potential. Their results are, how¬ 
ever, limited to shorter time scales than those obtained in kMC simulations. 
Both these methods must be applied for specihed potential or jump rates 
and at specihed temperature. By changing any of these, one has to perform 
the simulation again. Below we derive analytic expressions for dihusion coef- 
hcients for adatoms and dimers based on the existing picture of the potential 
energy landscape. Otherwise than MC or molecular dynamic calculations, an¬ 
alytic formulas express dihusion coefficient as a function of temperature and 
other model parameters [^22p23|l24] . It gives much more insight into the real 
physics of the system. Moreover, if new, more accurate data about the energy 
landscape are obtained, it is possible to improve the dihusion coefficient values 
on using the same formulas. Among theoretical approaches there are methods 
based on Langevin equation, Fokker-Planck equation or path integral formal¬ 
ism EH- Our approach, on the other hand, is built on the master equation. 
Similarly to the method developed by Titulaer and Deutch [25], which was 
later applied in [261127] for metal surfaces, we hnd the dihusive eigenvalue of 
the rate matrix. However, we do that by means of a variational method. 


In this work we calculate the dihusion coefficient for Cu monomer and Cu 
dimer on Cu(lll) and Ag(lll) surfaces. In such a way dynamics of homo- and 
heteroepitaxial system is compared for the same adatom. Such systems have 
been investigated both experimentally [HHini and theoretically [T^13II14II15II16II 
In addition to dihusion of a dimer, there have been also studies concerning 
adatom clusters of larger size [T^14lll5ll2nil28j . Hetero-dihusion appears to be 
easier with lower activation barrier and homo-dihusion is dominated by rela¬ 
tively high energy barriers between neighboring cells, whereas rotation inside 
each cell can run freely. 


We calculate dihusion coefficient using adsorption energies and energetic bar¬ 
riers for jumps between sites obtained in Refs. [T3][T6] by means of ab-initio cal¬ 
culations. We get analytical expressions for dihusion coefficients from which 
general relations between diherent dihusion modes at various temperatures 
are derived. Ehective dihusion coefficients are represented by prefactor and 
activation energy parameters. In general both these parameters show temper- 
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ature dependence which is in contrast with the global Arrhenius behaviour 
of the diffusion known for simple systems. Such temperature dependence of 
diffusion could explain some discrepancies in results obtained from different ex¬ 
periments [2^301131] ■ Analysis of diffusion coefficients also shows that whereas 
dimer diffusion on Ag surface happens via translation together with rotation, 
on Cu surface both moves are separated, and fast rotation is limited only to 
cells, whereas slow intercell jumps decide about rate of translational diffusion 
of dimers. 


2 Calculation of diffusion coefficients 


Variational approach to the calculation of the diffusion coefficient from a set of 
master equations was proposed in Ref. [32j. Since then it has been applied to 
many different systems of adatoms jumping on a crystal surface. So far mostly 
strictly theoretical systems have been considered. However, recently we have 
succesfully applied the variational approach to the diffusion of a single Ga 
adatom on the GaAs(OOl) surface c(4x4) a and /3 reconstructions |33] using 
the energy landscapes from Refs. [3l] and [35]. It is also possible to solve a 
problem of correlated particles using our approach. Till now single adatom 
has always been treated as a simple particle, that is a particle whose state 
could be completely described by giving its position on the lattice. It could 
jump to one of the neighbouring sites on the lattice provided it wasn’t already 
occupied. Such an adatom didn’t possess any additional degrees of freedom 
like for example orientation or shape. 

In this paper we describe dimer diffusion, namely Gu dimer on Gu(lll) and 
Ag(lll) surfaces. A dimer could be thought of as a compound particle or a 
cluster, that is an ensemble consisting of a certain number (two in the case of 
dimer) of simple particles that are bound together. Therefore, the move of the 
dimer should always be considered as the move of the whole, not as that of 
the individual constituents. The state of the dimer can be described fully by 
giving the position of its centre of mass, its length (the distance between the 
constituents) and its orientation along the underlying lattice. Therefore, the 
approach can be generalised to compound particles without any modifications. 

Diffusion coefficient of a single particle jumping on a crystal lattice is defined 
as the mean square displacement over time [T1I31I211124] 

Dnm = hm 4 < Ar„(t)Arm(t) > . (1) 

Indices n and m denote one of the space coordinates (x or y) on the lattice. 
Arn{m){t) is the displacement of the particle’s centre of mass in the n(m) di¬ 
rection after time t. Since we follow the centre of mass, it doesn’t matter if by 
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a particle we mean a simple point particle or a cluster. Position of its centre of 
mass is always well defined. The above definition is used for extracting the dif¬ 
fusion coefficient from the data obtained from scanning tunneling microscopy 
or molecular dynamics methods. 

In order to analyse diffusive motion on a crystal surface we can divide ar¬ 
bitrarily the lattice into unit cells which repeat themselves periodically over 
inhnite surface. Origin of each cell has position fj in the space and within 
given cell there are m adsorption sites, whose locations relative to the cell’s 
origin are given by aa, where a = 1,.. .m. Therefore, each adsorption site on 
the lattice is specihed unambiguously by parameters (j,a) and its position in 
the space is equal to = fj + da- Physically, position of a given adsorption 
site corresponds to the position of the centre of mass of the particle (simple 
or compound) occupying that site. A particle occupying the site (j,a) has the 
adsorption energy E{a). In general, it is possible for two or more different 
adsorption sites to have the same position. They will be distinguished by the 
index a though. For the coverage p (ratio of adparticles at the surface to the 
total number of adsorption sites) the probability of hnding a particle in the 
site of the energy is equal to Peq{a) = pexp[—(3E{a)]/J2'y^w[~PE{'y)]. 
The summation is done over all the sites in the cell and f3 = l/ksT is the 
inverse temperature factor. 

Diffusive motion on a surface consists of a series of thermally activated jumps 
between adsorption sites. According to the transition state theory, the jump 
rate (number of jumps per unit of time) is a function of the energy barrier 
along the path of the jump [36]. For a single jump from the (j,a) site to the 
(1,7) site it is equal to 

M^(j,a;^,7) = = i/exp{-/?[E(j,Q;;/,7) - E{a)]], (2) 

E{a) is the adsorption energy of the initial site of the jump and E{j,a-,l,'j) 
is the energy at the saddle point between the two sites involved in the jump. 
These are the energies that are found by means of ab-initio calculations [T6] . 
Variable u is the diffusion prefactor, whose exact value is a subject for a 
separate study, however, its value is very often assumed to be of order of 
10^^/s. It is sometimes also called the attempt frequency and its physical 
meaning is the number of attempts of the particle to jump out of its site per 
unit of time. The exponential part in ([^, on the other hand, is the probability 
of a successful attempt. Time evolution of a system is described by the Master 
equation 

^P(j,«;f) = - W^,,P{j,a;t)]. (3) 

are jump rates dehned by ([^ and P(j,a;t) is the probability that the 
site (j,a) is occupied. We are looking for the diffusive eigenvector of Eqs (|^. 
We will calculate this eigenvector by the variational method. Details of this 
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Fig. 1. Fcc(lll) lattice with two possible energy minima for a Cu monomer. 


method are described in Appendix A. 


3 Diffusion of Cu monomer and dimer on C'm(III) and Ag{lll) sur¬ 
faces 


Both copper and silver crystallize in the face centered cubic structure and their 
(111) lattices consist of two types of sites called fee and hep, and arranged 
in the way shown in Fig. The lattice constant for silver is equal to 4.09A 
compared to 3.6lA for copper. This difference in lattice constant is due to the 
interaction strength between substrate atoms and it also determines the shape 
of the energetic surface that is seen by adsorbed Cu atoms. 

The energy landscape for Cu monomer and Cu dimer on Cu(lll) was cal¬ 
culated in Ref. |T^ and on Ag(lll) in Ref. [H] using ab-initio methods. On 
the base of these calculations we can see that Ag surface observed by Cu 
adatom in hetero-diffusion process is more smooth. There is no evident differ¬ 
ence between jump rates inside and outside surface cell. At the same time in 
the homo-diffusion process Cu/Cu(lll) jumps from one cell to the neighbor¬ 
ing one are the slowest ones and determine the global diffusion rate. We use 
here the potential energy surfaces from Refs. [ms] to calculate the diffusion 
coefficient for both Cu monomer and dimer. For the dimer we calculate also 
rotational diffusion coefficient and compare it with the translational one. 

A single adatom can occupy either fee or hep site. These sites differ from each 
other by the position of the atoms lying directly below. An adatom occupying 
an fee site will therefore form bonds of different lengths and directions with 
underlying surface atoms than the one that occupies an hep site. It means 
that the potential energy of the adatom will depend on whether it resides at 
an fee or an hep site. However, according to both the experiment pain] and 
the calculations urns], the energy difference of these two sites is very small. 


5 











Fig. 2. Possible configurations of Cu dimer on the fcc(lll) surface. The sfh config¬ 
urations appear only on the Ag(lll). 

of order of single meV at both Cu(lll) and Ag(lll) surfaces. 

At the same time dimers formed at these surfaces see completely different 
landscape of potential energy. It is observed that dimers become quite stable 
and long-lived structures 09! and it is possible to study their diffusional 
movement. Dissociation of the Cu dimer on Cu(lll) has been never observed 
experimentally or in simulations. According to [I2| the energy barrier for such 
a process is 450 meV, which is much higher than for most of the possible jumps 
of the dimer. On the other hand, dissociation of the Cu dimer on Ag(lll) 
has occured at 700 K in molecular dynamics simulations [T6|, which suggests 
that the corresponding energy barrier is reduced on that surface compared to 
Cu(lll). However, the dissociation was only temporary and after a few ps the 
constituents recombined forming the dimer again. Therefore, the dimers can 
be treated as stable particles in our temperature range. 

According to Ref. [15] Cu dimer on Ag(lll) can be found in one of the four 
possible energetic minima (see Fig. |^. They differ by the type of site occupied 
by each single particle in the dimer (fee or hep) and by the distance between 
those particles. Different dimer conhgurations are called ff (two fee sites), hh 
(two hep sites), sfh (short fh) and Ifh (long fh). The ff and hh states differ 
only by the type of occupied sites while the distance between the particles 
is the same in both cases. On the other hand, in sfh and Ifh states, particles 
occupy the same pairs of sites but the distance between them can be either 
short (sfh) or long (Ifh). At Cu(lll) surface Cu dimer can be found in ff, hh or 
Ifli conhgurations only. The sfh state is not built on this lattice [T2|. Each of 
the dimer’s energy minima on this surface can appear in three variants which 
differ by their relative angular orientation. All the possible conhgurations of 
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Fig. 3. Diffusion coefficient for Cu monomer and Cu dimer on the Cu(lll) (top) 
and Ag(lll) (bottom) surfaces as the function of the temperature, v = 10^^/s. 

the Cu dimer on the fcc(lll) lattice are shown in Fig. 

It has been argned [16] whether the sfli and Ifh states are actnally potential 
minima on the Ag(lll) lattice or jnst metastable states. They he only a few 
meV below the calcnlated transition states to the ff and hh states and the 
potential energy landscape aronnd sfh and Ifh states is almost flat. However, 
anthors treat these conhgnrations as separate states and in onr calculations 
below we shall also consider the sfh and Ifh sites as energy minima. In other 
situation all barriers for jumps between states should be calculated from the 
beginning and then they can be used in our formula giving similar results as 
the one shown below. 


3.1 Monomer 


Structure of the energy landscape as seen by Cu monomer is the same for both 
Cu(lll) and Ag(lll) surfaces. Both the energetic minima and the barriers 
are found at the same positions. The only difference is the exact values of 
the energies. Therefore, we can derive a general formula for the monomer and 
then substitute specihc values for both Cu(lll) and the Ag(lll) surfaces in 
order to calculate the diffusion coefficients. 
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Fig. 4. Effective activation energy for the diffusive motion of Cu monomer and 
dimer, and for the rotational diffusive motion of dimer on the Cu(lll) (top) and 
Ag(lll) (bottom) surfaces as the function of the temperature. 


As it is apparent from the first-principles calculations [T21IT6] and from the 
experiments pimn] , Cu monomers are slightly more stable at fee sites than 
at hep sites. For Cu(lll) surface the difference between energies at those two 
sites is equal to 5 meV na and for the Ag(lll) lattice it is equal to 6 meV [T6] . 
The energetic barrier for the jump from an fee to an hep site on Cu is equal to 
0.041 eV and to 0.075 eV on Ag. The barriers for inverse processes are 0.036 
eV and 0.069 eV, respectively. We use these data in order to calculate the 
diffusion coefficient for a single Cu adatom on the Cu(lll) and the Ag(lll) 
surfaces. The variational parameters in (20) are equal to 0 because of the 
system symmetry. Therefore, we just have a sum of two terms, each related 
to one type of jump, which can be also expressed as 


® hF/cc—^-/icplF/jcp—> jee 

2 ^Vfec^hep T ^^hep—^fee 


(4) 


The parameter a in the above formula (and in all equations in this article) is 
the length of the lattice constant of the underlying surface, which is also equal 
to the bond length of Cu or Ag. It is related to the bulk lattice constant A by 
a = 4^ what gives acu = 2.55A for Cu and uaq = 2.89A for Ag surface. 

According to (|^, the diffusion of Cu monomer on Cu(lll) and Ag(lll) is given 
by the same formula as this for the one-dimensional diffusion of a single particle 











in the potential with two alternating minima. The temperature dependence 
of the monomer’s diffusion coefficients is shown in Fig. for Cu surface in top 
panel and for Ag surface in bottom panel, in both cases plotted as top lines. It 
can be seen that monomer diffusion on Ag surface is slightly slower than this 
on Cu lattice. Assuming a local Arrhenius behavior of the diffusion coefficient 
D = Do exp(—it is possible to calculate the effective activation energy 
Ea = For the monomer it is equal to 

fcc-^hcp^E 

hcp^fcc + fcc^^fcc—^hcp 

^Vfcc—>-hcp T ^^hcp—>-fcc 

It is obvious that a quantity defined in such a way depends on the tempera¬ 
ture. Ea for monomer on Cu lattice is shown in top panel of Fig. drawn with 
the middle line and Ea for monomer on Ag lattice it is shown with the top line 
in the bottom panel of Fig. As seen in both cases Ea does not change much, 
it starts at lower temperatures from the value equal to the height of the larger 
of the two barriers and at higher temperatures it goes down approaching the 
value which is the average of both barriers. Therefore, for Cu(lll) surface 
the effective activation energy goes from 41 meV to 39 meV. We can com¬ 
pare then with low temperature experimental value 37±5 [TU]. For Ag(lll) 
surface activation energy starts at 75 meV at low temperatures compared to 
65±9 meV in the experiment [8], then it goes down reaching 72 meV at high 
temperatures. 

Knowing the effective activation energy, it is also possible to calculate the 
effective prefactor Dq = D exp{f3Ea). Its behaviour is shown in Fig. [^with the 
bottom line in the top panel and with the top line in bottom panel. In both 
cases at low temperatures they decrease rapidly up to about 100 K and then 
remain constant at higher temperatures. 


3.2 Dimer diffusion 


For Cu dimer on the fcc(lll) surfaces the situation is more complex than for 
the monomer. As discussed above we have the following dimer configurations: 
ff, hh, fh that differ by the positions of two atoms of dimer. As we shall see 
below the distance between the atoms also matters and in general we will 
have sfh (short fh) and Ifh (long fh) positions. Moreover each configuration 
occurs in three different orientations. All the possible configurations found on 
the fcc(lll) lattices are shown in Fig. The map of the possible transitions 
between all those sites is quite complex and differs significantly between Cu 
and Ag surfaces, therefore, below we shall discuss both cases separately. The 
energy landscape and values of energy for both the sites and the barriers have 
been taken from Ref. [12] (for Cu(lll)) and from Ref. [IB] (for Ag(lll)) and 
used below to calculate the diffusion coefficients of Cu dimers on those lattices. 


9 





45 

40 

_35 

0 ) 

30 

^ 25 
o 

:=.20 

o 

“ 15 
10 
5 


Cu/Ag(111) 


Monomer 


kPimer 


0 100 200 300 400 500 600 

TTKl 


700 


Fig. 5. Effective prefactor for the diffusive motion of Cu monomer and dimer on the 
Cu(lll) (top) and Ag(lll) (bottom) surfaces as the function of the temperature. 
u = 10^^/s 


3.2.1 Cu dimer on Cu(lll) lattice 


Only three types of equilibrium states from Fig. are found for dimer on 
Cu(lll) lattice [TOlIT^ . These are ff, hh and Ifh. Conhguration sfh has too 
high energy at this lattice hence it cannot be used by dimer to realize intercell 
jumps. Energies of dimer states and barriers for transitions between those 
states has been worked out in Ref. na. We plotted the map of all the possible 
jumps in our system in Fig. Yellow circles represent adsorption sites for the 
dimer (its centre of mass on the lattice) and the lines represent transitions. 
Three Ifh sites with different subscripts have actually the same position but 
have been drawn as separate ones for clarity. The black lines represent intracell 
jumps, which are limited to a hnite area on the lattice. The dashed black line 
between ff and hh sites means that this jump is directly between these sites and 
does not pass through the Ifh site. Blue, red and green lines represent intercell 
motion in different directions, that is the one which goes beyond given cell. 
Energy landscape along the path marked in red is plotted in the lower part 
of Fig. Plotted path leads dimer from one cell to the next one. Each cell 
is visible in top panel of Fig in the form of star. Intercell jumps have much 
higher energy barrier than intracell ones. In our calculations we shall take 
into account all the intracell jumps and the lowest in energy intercell motion, 
namely the concerted sliding between ff and hh sites, which has an energetic 
barrier of 120 meV. It is shown in Fig. [^how this particular intercell jump is 
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Fig. 6. Scheme of the transitions on the Cu(lll) (top) surface and energy profile 
(bottom) along the marked path for a single Cu dimer. Names of the sites are 
assigned according to Fig. 

realized and how it differs from the intracell jnmp of the same type. Barriers 
for all the remaining intercell jnmps (inclnding the dimer’s dissociation) are 
at least three times larger, therefore these jnmps are highly nnlikely and will 
be omitted in fnrther considerations. 

We denote rate of the intercell jump by Vff^hh, while the intracell jumps 
will have rates Wfj^hh, Wff^ifh, Whh^ifh and Wifh^ifh- Symbols with tilde 
over them will be equal to corresponding rates multiplied by the equilibrium 
probability of the initial site’s occupancy, for example Vff^hh = PilVfS —^hh * 
We shall use this notation in the whole article. 

It can be seen in Fig. |^that the system possesses three axes of symmetry (one 
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Fig. 7. Realizations of intercell jumps via hh—7-ff jumps at Cu(lll) surface (top) 
and via sfh state at Ag(lll) surface (bottom). 


Cu/Cu(lll) path 2 
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Cu/Cu(lll) path 1 
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Fig. 8. Different diffusion modes on Cu(lll) (top) and Ag(lll) (bottom) surfaces. 
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Fig. 9. Rotational diffusion coefficient for Cu dimer on the Cu(lll) (top) and 
Ag(lll) (bottom) surfaces as the function of the temperature, v = 10^^/s 



Fig. 10. Effective prefactor for the rotational diffusive motion of dimer on the 
Cu(lll) (top) and Ag(lll) (bottom) lattices as the function of the temperature. 
V = 10^^/s. 

along the a sites, one along the b sites and one along the c sites). The lattice 
sites labeled by different subscripts (a, b or c) are energetically equivalent, that 
is they have exactly the same potential energy and the same neighboring sites 
with the same jump rates to them. In order to apply our variational method, 
we need to introduce some parameters that depend on the state of the system. 
The state of the system is given completely by specifying in which of the sites 
shown in Fig. |^the dimer currently resides. Each of those sites should have 
its own geometric vector phase, which can be written for convenience in the 
angular form 6 = dr (cos sin d^). The form of these phases should reflect the 
symmetry of the potential felt by the dimer at given site. For site a it can be 
seen that there is symmetry with respect to the y axix, therefore we should 
have d^ = |, which gives = (0, 5“). The angular parts of the corresponding 
b and c phases should then be equal to ^ and respectively. Therefore, the 
only variational parameters will be the radial parts, that is d^-^, d^^ and d\.^^. 
From now on we will omit the subscript r. 
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Fig. 11. Scheme of the transitions on the Ag(lll) (top) surface and energy prohle 
(bottom) along the marked path for a single Cu dimer. Names of the sites are 
assigned according to Fig. [2j 


We insert such defined phases into the variational formula (20) and, using 
simple differentiation, minimise the expression with respect to the parameters. 
Resulting expression for phases are given in Appendix B. Finally the diffusion 
coefficient for Cu dimer at Cu(lll) surface is 


D --aWff^hh 


{l + 3 Wff^hh{Wff ^Ifh + Whh^lfh) 


+2Wff^ifh^hh^ifh\ {Wff^hh + 2V//_j,/i/i) 
+ Whh^lfh) + 2Wff^ifhWhh^lfh 


-1 


( 6 ) 


As expected, the intercell jump rate Vff^hh is fundamental for the diffusive 
motion of the dimer. The hrst term in the above sum describes the diffusion 
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controlled exclusively by the // —)■ hh intercell jumps, which are quite slow. 
However, the second term has a large contribution to the total value of the 
diffusion coefficient. It describes the motion which is partially realised by the 
much faster intracell jumps. Therefore, we see that even though those jumps 
are unable to trigger the diffusion by themselves, they can enhance it signif¬ 
icantly by bypassing some of the slow intercell jumps. These two modes of 
diffusion: one omitting intracell jumps given by hrst term of sum ([^ and the 
second through the inter and intracell jumps are illustrated in top panel of 
Fig. 1^ Fig. shows diffusion coefficients for both modes. It is seen that the 
one plotted with the bottom line is much slower than the second one what 
means that appropriate diffusion path is less probable. 


It should be also noted that jump rate Wifh^ifh is not included in the above 
expression. It is understandable, since those jumps only rotate dimer between 
its different angular orientations. They do not move the dimer centre of mass at 
all. However, it is also possible to calculate the rotational diffusion coefficient. 
By analogy to Q it is dehned as 

D, = lim - < 02 > . (7) 


The variational formula (20) can be used in the same form, if instead of the 
translational degree of freedom f we use the angular one 6. The variational 
phases coupled to the angle 9 go to zero, and it is easy to understand as all the 
sites are symmetric with respect to change of the angle. Thus the rotational 
diffusion coefficient for the Cu dimer on the Cu(lll) is equal to 


Dr — yy + ^hh^ifh + ‘2Wifh^ifh)- (8) 

io 

It implies that the rotational diffusion is realised in three independent ways. 
The contribution of the Ifh —)■ Ifh jump is twice as large as these of // —)■ Ifh 
and hh —)■ Ifh because there is only one Ifh —)■ Ifh jump needed in order to 
rotate the dimer by 120 degrees. Both // —)■ Ifh and hh —)■ Ifh need to be 
performed twice for the same effect. 


It can be seen in Fig. that the dimer diffusion coefficient is much smaller 
than that for the monomer. It is more difficult to compare translational and 
rotational motion of dimer because one is measured in nmf /s and the second 
in radf/s. In order to compare them both we multiplied Dr by (P/4:, square 
of half of the dimer length. Such value can be attributed to the diffusion of 
one of the atoms of dimer in the rotational movement. In Fig. the upper 
line shows Cu dimer rotation in this case and it can be seen that it is as fast 
as monomer diffusion over the surface and at the same time much faster than 
dimer translational diffusion. It means that dimer rotates at given cell and 
this rotation does not couple to the translational diffusion. We also calculated 
the activation energies for these motions just as we have already done for the 
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monomer. The formula for the translational dimer’s motion is too complex to 
write it down, therefore we just plot the result in Fig. The expression for 
activation energy for the dimer’s rotation is 


Ea = {Eff + AEff^ifh)Wff^ifh + {Ehh + AEhh^ifh) 
Whh^lfh + ‘^{Eifh + AEifh^ifh)Wifh^ifh / 

+ ^hh^lfh + ‘^Wifh^ifh)— < E > 


(9) 


where 


<E>= 


^ lA ^ — 

Eff exp{-/3Eff) + Ehhexp{-/3Ehh) + Eifhexp{-(3Eifh) 
exp{-l3Eff) + exp{-l3Ehh) + exp{-^Eifh) 


( 10 ) 


is mean value of dimer energy in the equilibrium conditions. The hrst term 
in is the weighted arithmetic mean of the energies at the saddle points over 
which the dimer jumps while changing its orientation. The second term is the 
average energy that the dimer has at a given temperature. We show all acti¬ 
vation energies for Cu(lll) surface in the top panel of Fig. We see that the 
activation energy for the dimer translational motion is the largest one and this 
for dimer rotation is lower even than this of the monomer diffusion. Therefore 
dimer motion consists mostly of rotation, while the probability for the move 
of the center of mass is much lower. Rotation is realized via intracell jumps 
and we can compare the calculated value Ea=20 meV with measured at low 
temperatures Ea intracell diffusion as 18±3 meV [10]. Monomer diffusion is 
somewhere in between. We can also see how prefactors for translational dimer 
motion (Fig. and rotational motion (Fig. 10) depend on temperature. It can 
be seen that prefactor for dimer diffusion is higher than that for monomer, but 
decreases with temperature in the same way, while prefactor for rotational dif¬ 
fusion increases with temperature reaching constant value from below. Below 
we will compare this behavior with this on Ag surface. 


3.2.2 Cu dimer on Ag(lll) lattice 


For Cu dimer on Ag(lll) there is one additional lattice position, namely the 
sfh site (short fh), which wasn’t present in the homoepitaxial case. According 
to Ref. [16] there are four types of transitions between Cu/Ag(lll) sites. 
The geometry of jumps on this lattice is signihcantly different than that on 
Cu/Cu(lll) and it requires a separate consideration. The map of the energy 
minima and the possible jumps in the system is shown in Fig. [T^ The path 


along which we plot energy potential is shown in top panel of Fig. H It can 
be seen that it leads from one cell to the next one. As before three neighboring 


cells can be seen in Fig 11 Three Ifh sites in each cell have the same position of 
the center of mass despite being drawn separately and the dashed line between 
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Fig. 12. Logarithmic scale plot of diffusion coefficient as a function of temperature 
(left) and as a function of inverse temperature (right) for Cu monomers and dimers 
on the Cu(lll) and Ag(lll) surfaces at low temperatures. Curves are described 
from top down. Two lines in the middle show monomer and dimer at Ag surface 
V = 10^^/s. 

// and hh sites represents the direct transition between these sites. 

We should point out some differences between the homoepitaxial and the 
heteroepitaxial case. First of all, as we have already mentioned, sfh position 
is added. It means that the direct intercell // —)■ hh jump is replaced by jumps 
through the sfh site (Fig. [^. Because of this there exist more sequences of 
jumps that are possible, since the dimer from the sfh site can jump either 
to one of the two // sites or to one of the hh sites. Therefore, it is difficult 
to state unambiguously which jumps are intracell and which are intercell. 
Another difference is lack of hh ^ Ifh and Ifh —)■ Ifh connections. 




We apply the variational method to the hetoroepitaxial system in the same 
way as we did for the homoepitaxial one. First we assign geometrical phases 
to each of the sites taking into account the symmetry, which is the same as 
in the previous case. Then we insert phases into Eq. (20) and minimize this 
expression with respect to them. Resulting formulas for phases can be found 
in Appendix B. And hnal expression for diffusion coefficient is 


D — 3a‘^{Wff^hhWff^sfhWhh^sfh{Wff^sfh + Whh^sfh) 

+ 2^ff^sfhWhh^sfh{Wff^hhWff^ifh 

+ ^^f f^lfhWhh^sfh 

+ Wff^sfhWhh^sfh) + 

+ Wlj^^gfff)}/{2Wff^sfhWhh^sfh{Wff^hh + Wff^ifh 

+ ^ff^sfh + Whh^sfh) + 

+ + Wff^ifhWff^sfh + Wff^ifhWhh^sfh) 

+ ( 11 ) 

The above formula is quite complex, however, there are still a few things that 
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can be deduced from it. First of all, we can see that unlike the Cu(lll) surface, 
there is no single type of jump on Ag(lll) that is crucial for the diffusion. 
On the other hand, the diffusion needs at least two different types of jumps in 
order to occur. One single type of jump is unable to trigger diffusion by itself. 
However, it seems that the // —)■ sfh jump is the most important one since 
it appears in almost all the diffusion paths. As we can see, the distinction 
between intracell and intercell jumps is not so obvious here as it was on the 
Cu(lll) surface. 


When we ignore the two slowest jumps i.e. Wff^hh = Whh^sfh = 0 the diffu¬ 
sion coefficient formula simplihes to 


4 Wff^sfh + Wff^ifh 

It corresponds to the zigzag motion observed by STM in Ref. [8]. 


( 12 ) 


For a more complex case, where the least probable transition is still Wff^hh = 
0 but Whh^sfh 7^ 0, the diffusion coefficient is 


D = 3a^ 


IF. 




{w 


114 


ff^lfhyVhh^sfh 


+ IF 


Wf 


ff^sfhyy ff^ifh 


3W 


ff^lfhWhh-^sfh 


f f^sfhWhh^sfhj / 

(iT/z^s/h + Wff^ifh + ^hh^sfh) 


(13) 


Here, in addition to the zigzag motion, we have also two other ones. It is also 
seen that jumps and Whh^sfh by themselves are not able to make 

the dimer diffuse along the surface. It is understandable because these jumps 
are between two completely different pairs of states. In the bottom panel of 
Fig. [^we plotted both approximations Eqs. ( I2p!3 ). In Fig. |^the successive 
dimer paths are shown. It can be seen that it is not zigzag motion (plotted in 
the lowest dashed line) alone that is responsible for the dimer diffusion. With 
two more modes added to this simplest one we can reproduce diffusion curve 
more precisely. The most probable dimer movement can be quite complicated 
as shown in Figj^ 


The effective activation energy in the most general case is expressed by a 
complicated formula, therefore we do not write the equation, but just plot the 
results in Fig. Its behavior is a bit more complex than for the monomer’s 
case, what is not surprising since the dimer’s map of transitions is also more 
complex as is the diffusion coefficient. At low temperatures the activation 
energy for the dimer is 75 meV compared to 73 meV in the experiment |H] 
and the molecular dynamics calculations [16]. Between 0 and 100 K activation 
energy for dimer diffusion is higher than this for monomers and for higher tem¬ 
peratures Ea for dimers decreases below value of Ea for monomers. Effective 
prefactor for the dimer plotted in Fig. is lower than this for monomers what 
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compensates the difference in activation energy. At lower temperatures where 
Ea for dimers is higher than that for monomers both diffusion coefficients have 
the same value as seen in Fig. 12 This is different than at Cu surface where 
even at low temperatures dimer diffusion coefficient is much lower than this for 
monomers. It can be also seen that at low temperatures dimers at Ag surface 
move much faster than at Cu surface. The same diffusion coefficient is shown 
in the right panel of Fig. 12 In this scale all lines are straight. Each curve has 
different slope, given by activation energy. However, it is impossible to see any 
difference in activation energies as a function of temperature in such plot. 


The rotational diffusion coefficient for the Cu dimer on Ag(lll) lattice is equal 
to 

Dr = + Wff^sfh + Whh^sfh) (14) 

It implies that the rotational diffusive motion in our case is a sum of three 
parallel motions, each of them related to a different type of jump. Transition 
Wff ^hh is not present in the expression, which is understandable because 
this transition does not contribute to the change of the dimer’s orientation 
as do the other three transitions. In contrast to the translational diffusion, 
the rotational one requires just one type of jump in order to occur. It is with 
agreement with a simple analysis of the dimer’s angular motion on the lattice. 
Each of the jumps whose rates appear in the expression can by itself rotate the 
dimer by the full angle. As was in the case of Cu/Cu(lll) it is not possible to 
compare diffusion coefficients for different types of motion (translational and 
rotational) directly, however, as in previous case in Fig. we plotted value Dr 
multiplied by square of the half of dimer length. It can be seen that rotation 
on Ag surface is much slower than rotation on Cu and it is very close to the 
translational diffusion of dimers on Ag. It can be understood that particles 
rotate while they move forward. 


As was in the case of the translational diffusion, we can also calculate the 
effective activation energy for the rotational diffusion 

Da = [{Eff + AEff^ifhWff^ifh + {Eff + AEff^sfh) 

^ff^sfh + {Ehh + AEhh^sfh)Whh^sfh\ 

+ Wff^sfh + Whh^sfh]- < E > (15) 

where < E > is mean energy of the dimer in the equilibrium state averaged 
over all four possible states. Note that even if not all transitions on the lattice 
are engaged in the rotational diffusion, the energy average is taken over all 
possible energies in the system. 


Because there are three parallel paths of rotational diffusion, the dimer can 
choose any of them independently. At low temperatures the effective activation 
energy for the rotational diffusion is 72 meV, which is also the energy barrier 
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for the // —)■ sfh transition. Even though the barrier for the hh —?■ sfh 
transition is even lower, this transition is less favourable because the hh site 
lies 13 meV higher than the // site and consequently, the saddle point for 
the hh —)■ sfh jump lies higher than for the // —)■ sfh jump. Therefore, the 
// —)■ sfh jumps make the easiest diffusion path. The effective activation 
energy then slightly rises reaching 73 meV at 50 K, possibly due to the other 
jumps becoming more important. Then above that temperature, as the dimer 
gets more and more energy, the activation energy starts to fall down very 
quickly reaching 


ET = 


‘^Eff T AT/ 




+ AE 


ff^sfh 


+ Ehh + AE, 


hh—^sfh 


Af 


Ehh + Eg 


« XmeV 


(16) 


in the limit of inhnite temperature. 


In Fig. IT we show an effective prefactor for the rotational diffusion. Its be¬ 
haviour is completely different than that of the translational prefactor. First it 
rises very quickly and about 50 K reaches a maximum slightly above 6T0^^/s, 
which is almost twice as high as the value at 0 K. After that maximum the 
prefactor goes down, just as in the translational cases. 


The different behaviour of both the activation energy and the prefactor for 
the rotational motion compared to the translational one can be explained due 
to the fact that the former is more parallel while the latter (especially for the 
monomer) is more serial. 


4 Conclusions 


Systems of diffusing Cu adatoms and dimers were compared at (111) surfaces 
for homo and heteroepitaxial systems. Cu and Ag crystals have different lat¬ 
tice constants and as a result potential energy surface seen by adatoms varies. 
However, for both surfaces diffusion of single adatoms is simple and has the 
same character and values given by appropriate activation energies for these 
processes. For dimers the situation is completely different. At Ag surface even 
if the lattice of dimer jumps is quite complicated, involving all types of possi¬ 
ble moves, all jumps are of the same rate and as an effect at low temperatures 
diffusion is quite fast. At higher temperatures it slows down. For homogeneous 
diffusion of Cu dimers the energy surface is more complex. Jumps inside hexag¬ 
onal cells are much faster than these between cells. Diffusion either happens 
only via intercell jumps - this is slower path, or it partly goes via intracell 
jumps, which is slightly faster. Both diffusion channels together give diffusion 
coefficient that is lower than this for Ag surface, but at higher temperatures. 
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around 600 K the relation turns back. It is interesting to compare tempera¬ 
ture dependence of rotational movement in both cases. They both have rather 
simple form being a sum of independent diffusional modes, but their effective 
activation energies are very different. Whereas activation energy for rotational 
motion at Cu surface is very low and almost independent of temperature, Ea 
at Ag is higher and decreases with temperature. As an effect rotation at Cu 
surface is as fast as single particle movement at this surface and it is localized 
inside cells, whereas rotation at Ag lattice is as fast as dimer movement and it 
happens simultaneously with translational movement of dimers. At low tem¬ 
peratures it is the same as monomer diffusion at Ag surface. It can be seen 
that even if both surfaces have similar geometry and not very distinct lattice 
constants, the character of resulting monomer and dimer diffusion is different. 

The general consequence of diffusion over complex energy surface is tempera¬ 
ture dependence of activation energy and prefactor of diffusion coefficient. For 
such an effect to appear there is no need for special interactions or other cor¬ 
relations in the system. Interestingly prefactors in our case drastically change 
on coming from low to high temperatures as shown in Fig 5. Such tempera¬ 
ture dependence of prefactors for monomer or dimer diffusion could explain 
some discrepancies in results obtained from different experiments [2^301131] . 
Especially change of prefactors ratio for diffusion at step and on terraces for 
Ag/Ag(lll) system [2^11^ can be considered as an effect of multistate par¬ 
ticle diffusion. We have shown above that depending on the temperature the 
character of diffusion on smooth (111) metallic surfaces changes due to small 
differences of the energy barriers for adatom or dimer jumps. We have also 
shown that at these surfaces dimer diffusion is more complicated, so it can be 
regarded as a good candidate for the interpretation of the experimental data. 
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6 Appendix A 


Because we deal with an inhnite periodic lattice it is convenient to take the 
Fourier transform of the equation ([^ 

^.Pa{k] i) = ^( 7 ; (y)P^{k] t) + M{a; a)Pa{k; f), (17) 
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where Pa{k;t) = exp(iA;f^)P(j, a; t). Elements M( 7 ; a) are expressed in 
terms of the jump rates in the following way: 

M{a-,a) = 

M{r, a) = ( 18 ) 

Sum in the hrst expression is over all possible transitions from the state a to 
the neighboring states. The second expression is written for each transition 
from 7 to a and I' = I as long as jump is within chosen cell, whereas /' = / ± 1 
when jump of adparticle transforms it to the neighboring cell. 

The matrix M contains all the information about system dynamics. Its eigen¬ 
values correspond to dynamic modes that describe relaxation of the system 
towards the equilibrium. It can be shown that of all the eigenvalues the diffu¬ 
sive one has the smallest absolute value pni25ll32j . Therefore, we can use the 
variational method in order to obtain the diffusion coefficient 

wMv ^ 

lim Xd = lim —= —kDk, (19) 

|fc|-5-0 |fc|-!>0 WV 

where w and v are left and right trial eigenvectors, respectively. Their elements 
are related to each other by w*^Peq{oi) = Va- Assuming Wa = where 0 q, 

is the geometric phase of the adsorption site a, we get the variational formula 
for the diffusion coefficient in the form 

kD^ark = lim ^ 

|fc|—^0 Q>7 

= ^ W^,,P,q{a)[k{rl + (20) 

Q!>7 

The above choice of variational vector has been so far valid in the cases of 
single adatom diffusion. Each of the vector phases corresponds to exactly one 
adsorption site within a single cell and each of its vector components is coupled 
either to x or y direction on the surface. The same approach can be used to 
describe the diffusion of the dimer center of mass. Moreover, the procedure can 
be also applied for rotational degrees of freedom instead of the translational 
one in order to calculate rotational diffusion coefficient. 

The total number of variational parameters, whose values should be chosen 
in such a way that the above expression is minimized, is twice the number of 
different types of adsorption sites. However, since only the phase differences 
contribute to the expression, we can always shift all the phases simultaneously 
by the same value. In such a way we can set one of the vector phases to 0. 
Moreover, the remaining phases can be related to each other due to the sys¬ 
tem’s symmetry, which leads to further reduction of independent parameters. 
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Ultimately we have a set of several linear equation, one for each independent 
parameter. We insert solution of these equations to Eq. (20) and receive hnal 
expression for the diffusion coefficient. 


7 Appendix B 


First step in the calculations of diffusion coefficient described above is to hnd 
variational parameters from the equation (20). Phases are associated with sites 
of the lattice that represents possible jumps in the phase space. Dimer con- 
hgurations are dehned by positions of both particles. In the case of Cu dimer 
diffusion at Cu(lll) surface sites are marked by //, Ifh and hh. Each site has 
attributed phase vector of orientation consistent with the lattice symmetry. 
Final expression for diffusion coefficient given by ([^ is obtained for Sifh phase 
which turns out to be zero and the other two 


Sff = 


w 


ff^hh 


+ Vt 


ff^hh) 


+ w. 


.W-, 


ff^hh 


hh^lfh\ 


2V, 


ff^hh) 


+ ^ff^lfhWhh^lfhiWff^hh — Vff^hh 
+ + ^hh^lfh + ‘^^Ifh^lfh) 

+ Wff^ifhWifh^ifhiWff^hh + 2Vff^hh) 

+ ^hh^lfhWlfh^lfhiWff^hh “ 

+ 2 Vff^hh){Wff^ifh + Whh^lfh) 

+ ‘^Wff^ifhWhh^ifh] [Wff -^Ifh + Whh^lfh 
+ 2Wifh^ifh]} 


( 21 ) 


Ohh - {Vv - - - ZVff^hh 

+ Vff^hh) + Wff^ifhWhh^lfhiWff^hh — Vff^hh 

+ + Whh^lfh + ‘^Wifh^ifh) 

+ ^ff^lfhWlfh^lfh{Wff^hh - 4:Vff^hh) 

+ Whh^lfhWifh^lfhiWff^hh + 

/{[{^ff^hh + 2Vff^hh)(Wff^ifh + Whh^ifh) 

+ 2Wff^ifhWhh^ifh] [Wff^ifh + Whh^ifh + ‘2'Wifh^ifh]} ( 22 ) 


) + ^hh^ifh' 


■ lU. 


ff^hh 


When Cu dimer diffuses on Ag(lll) surface it has four different positions 
//, sfh, Ifh or hh. In the most general case, with all the jumps taken into ac¬ 
count diffusion is given by Eq. © where the following variational parameters 
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were used 


- -{^^ff^lfhWhh^gfh + 2^ff^hhW^h^^fh 

- Wff^^ff^Whh^sfh + Wff^sfhWff^ifhWhh^sfh 
1^2' 

- + Wff^sfhWff^hhWhh^sfh 

+ 2^ff^ifh{Wff^hhWhh^sfh + Wff^sfhWff^hh)} 

/ {‘^^ff^sfhW^h^^fh + 2^ff^lfhWlh^sfh 

+ + 2Wff^^^ff^Whh^sfh 

+ Wff^^fi^Wff^hh + 2Wff^sfhWff^ifhWhh^sfh 
+ 2Wf f^sfhWf f^hhWhh^sfh + Wf f^ifhWf f^hhWhh^sfh 

+ ^ff^sfhWff^ifhWff^hh} (23) 


3 1 

5hh = 

+ + 2^ff^sfhWff^ifhWhh^sfh 

3 ~ 2 ' 

- 2 ^ff^sfh^ff^hh - Wff^sfhWff^hhWhh^sfh 

~ 2^ff^ifh{Wff^hhWhh^sfh - Wff^sfhWff^hh)} 

/ {‘^^ff^sfhW^h^^fh + 2^ff^lfhWlh^sfh 

+ + ‘^W^ff^^fJ^Whh^sfh 

+ Wff^^fi^Wff^hh + 2Wff^sfhWff^ifhWhh^sfh 

+ 2Wf f^sfhWf f^hhWhh^sfh + Wff^ifhWff^hhWhh^sfh 

+ (24) 




Whh^sfh 


(1 + '^// ~ ~ 2(5fe/i — 


1 ^ff 
4 2 


1 

2 


(25) 
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For the case Wff^hh = Whh^sfh = 0 the phases simplify to 


hf - 


fF: 


ff^sfh 


- W 




Shh = 0 

^sfh = ^Ifh = 


2{Wff^sfh + Wff^ifh) 


2{Wff 

-^sfh 


(26) 


and we obtain Eq. (12). For a more complex case, where the least probable 


transition is still Wff^hh = 0 bnt Whh^sfh 7 ^ 0 , we have the following valnes 
of phases 


= 


ff^lfhWhh^sfh + 


-w 


Wf 


ff^sfhyy ff^ifh 


/ 


-Wff^ifhWhh^sfh 


+2Wff^sfhiWhh^sfh + Wff^sfh + Wff^ifh) 


(27) 


Shh — — 


rlF 


ff^lfhWhh^sfh 


+ 1F 


f f—^sfh^^hh—^sfh 


+ 2^ff^sfhWff^ifh 


/ 


-Wff^ifhWhh^sfh 


+2Wff^sfh{Whh^sfh + Wff^sfh + 


(28) 


Ssfh - 

Sifh = 


~ 2 ~ 

4 2 


(29) 


and these phases lead to Eq.(13) 
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